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The QCD phase transition is studied on 32 x 8 quenched lattices. One expects from simple arguments that 
the quenched chiral condensate should diverge as the chiral limit is approached. Previous studies have not been 
able to observe this behavior, perhaps due to relatively large lattice spacings used. The topic is revisited here on 
Nt = 8 lattices, where the physical spacing is approximately (2 GeV)^ 1 . 



1. INTRODUCTION 

Quenched QCD provides an interesting oppor- 
tunity to study the role of topological effects on 
the axial U(1)a symmetry above the chiral sym- 
metry restoring phase transition. Because it will 
turn out that our results elucidate features of the 
full dynamical theory as well, they have bearing 
on a wide class of problems, in particular the or- 
der of the QCD phase transition itself Q. 

We will begin by reviewing the role of small 
eigenvalues of the Dirac operator in producing 
a chiral condensate. We will then show that in 
the quenched approximation, where the fermion 
determinant is removed from the path integral, 
there is a possibility for a dramatic effect which 
would be seen at temperatures just above the de- 
confining transition. In this case, U(1)a break- 
ing effects manifest themselves as a divergence in 
the condensate because they are not regulated by 
suppression from fermion loops. This situation is 
particularly easy to understand in an instanton 
gas picture, so we will explain it in some detail 
there. 

We conclude by presenting results from our own 
investigations of 32 3 x 8 lattices on Columbia Uni- 
versity's QCDSP parallel supercomputer, and an- 
alyzing them in this context. 
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2. THEORY 

It is not difficult to show that the chiral conden- 
sate is just the trace of the the Euclidean fermion 
propagator (i P + m) . By looking at this trace in 
the representation of the eigenvectors of the Dirac 
operator itself, we can separate the contributions 
to the chiral condensate into two parts M: 
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Here the first term represents the contribu- 
tion of exact zero modes, iV zm in number, which 
are topological in nature while the second results 
from the continuum of Dirac eigenvalues. 

2.1. Instantons 

A simple picture which will make the role of 
these zero modes clear is one of dilute instantons. 
Though we will see below that more general con- 
siderations may apply, it is particularly instruc- 
tive to think for a moment in terms of this simple 
concrete model. 

To first approximation, we can take dilute 
to mean that there are no interactions of any 
kind between widely separated instantons or anti- 
instantons. In this case, the number of zero 
modes should go like iVj (the number of instan- 
tons) plus N a (the number of instantons). In such 
a picture, as the volume grows, so also does the 
number of zero modes, and the density of exact 
zero modes does not vanish in the infinite volume 
limit. 

We can extend these ideas to include a more 
plausible situation in which there are weak in- 
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teractions between the (anti-)instantons. Even 
though the presence of interactions will give a 
topological charge v — Ni minus N a , and thus 
only exact zero modes, there will still be nearly 
zero modes A nZ m shifted away from A = only 
slightly by the presence of the interactions. So, 
although the true zero modes in this case will be 
suppressed like ^= in the infinite volume limit, 
the nearly zero modes will persist. 

2.2. Expectations 

We might not expect such configurations to 
contribute significantly to an observable such as 
the chiral condensate. Even though the contribu- 
tion appears to be divergent in the chiral limit, 
the fermion determinant det(i p+m) is suppress- 
ing the contribution of such configurations with 
a factor going like m zm . In a quenched simula- 
tion however, we effectively set det(i p + m) = 1. 
This will eliminate the suppression of such con- 
figurations, amplifying their contribution to the 
chiral condensate, and thus manifesting the x j m 
divergence which one would naively expect from 
(!)■ 

The instanton gas model leads us to expect a 
finite density of instantons in enough configura- 
tions to have weight in the path integral. An 
immediate consequence is that the quenched ap- 
proximation will result in a I /m divergence in the 
chiral condensate. Even in the presence of weak 
interactions between the instantons, all of this 
will continue to hold true, so long as the masses 
being studied are greater than A nzm . 

We can see however, that the situation is based 
on a more general idea than the dilute instan- 
ton model alone. Though such a picture gives us 
a concrete way of thinking about the situation, 
it is possible that the x j m behavior seen in that 
specific approximation is really telling us some- 
thing about the general structure of the path in- 
tegral itself. Namely, the lack of suppression of 
the fermion zero modes by the quark determinant 
may be itself sufficient to produce this divergence. 

3. DESCRIPTION OF CALCULATION 

Previous quenched calculations in this temper- 
ature region by the Columbia group || , found that 



chiral symmetry remains broken when the config- 
urations take on either of the complex phases of 
the Wilson line. Because we are ultimately inter- 
ested in an investigation relevant to unquenched 
physics, we restrict our attention to the purely 
real phase, in which chiral symmetry is restored 
above the deconfining transition. 

If we hope to resolve topological effects such as 
instantons, we must insure that the lattice spac- 
ing is sufficiently fine that these features are accu- 
rately represented. If the features we are looking 
at are indeed instantons, we would expect that 
the small ones would begin to decouple as a result 
of asymptotic freedom. Taking the rho Comp- 
ton wavelength as an approximate intrinsic scale 
of QCD, it is objects of this size that we might 
expect to need to be able to resolve. N t =A for 
ft ~ ftc gives only 

^lattice ~ 

(I Ge^)- 1 . Cutting 
the spacing in half by going to N t = 8, gives us a 
lattice spacing closer to (2 GeV)^ 1 , and increas- 
ing confidence that we are probing the interesting 
physical region. 

In a quenched simulation, we are free to mea- 
sure a variety of valance masses on the same en- 
semble of configurations because the mass plays 
no role in the generation of the ensemble. We do 
not however choose to measure extremely light 
masses, because the lattice itself provides a cut- 
off for the smallest available eigenvalues for i p. 
Masses below this cutoff are in effect all identi- 
cal for the purposes of computing the chiral con- 
densate. For confined physics, we know that the 
smallest eigenvalues have size set by the inverse 
total volume. This is the basis of the appearance 
of nonzero eigenvalue density around the origin 
and nonzero chiral condensate. This represents a 
severe compression due to QCD dynamics relative 
to the free Dirac eigenvalues, for which A m i n is of 
order ? where L is just the extent of the phys- 
ical volume along its longest side. This density 
decreases sharply as we increase ft through the 
transition, and continues to be slowly suppressed 
as we continue to yet higher temperatures. For 
this reason, y is a lower bound on A m i n . For a 
volume of 32 3 x 8 = 2 18 , we expect masses less 
than or equal to about m q — I0~ 6 to be uninter- 
esting for this reason. 

The -/V t =4 Columbia study mentioned above 



was performed at ft = 5.710. It is clear that 
running at unnecessarily high temperatures will 
suppress the signal we are looking to observe, 
while running too close to the transition intro- 
duces the new difficulty of tunnelings into the 
confined phase. For Nt — 8, taking the critical 
coupling to be Pc Nt=8 — 6-06 ||, the correspond- 
ing coupling would be (3 CN _ a + A/3 = 6.08. We 
chose to run here also to make comparison with 
this existing N t — 8 data as meaningful as possi- 
ble. 

3.1. Simulation Details 

This simulation was run on a 512-node pro- 
totype of the Columbia University 8192-node 
QCDSP supercomputer ||. The prototype sus- 
tains a peak performance of 25 Gigaflops. The 
evolutions were generated by a simple heat bath 
algorithm for pure gauge SU(3) (7]|| . The spec- 
trum data represents 50 measurements, separated 
by 250 sweeps of the pure gauge evolution, using a 
conjugate-gradient inverter to calculate the chiral 
condensate with a staggered fermion Dirac oper- 
ator. The source vectors for the CG inverter are 
generated randomly to produce a noisy estimator 
for the fermion propagator trace. All masses are 
measured using the same fermion source vector. 
This calculation was done on a 32 3 x 8 volume at 
[3 = 6.08. 
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N.=4 and N.=8 Data 
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Figure 1. For N t = 4 and Nt = 8, we see 
a power law fits the available data over a wide 
range of masses. 



gin. On the other hand, it is quite possible, that 
the opinion that there is validity to the prediction 
of this divergence beyond the dilute instanton ap- 
proximation, is itself severely misplaced. 



4. RESULTS 

We observe no signal for the expected %„ di- 
vergence. In fact for the mass range 10 -5 < m < 
5 x 10 -2 the curve is fit very nicely by a power 
law with an exponent of 0.91 (x 2 /dof = 2.0) . 
This can be compared to the Nt — 4 data, where 
a similar fit was found to well accommodate the 
data, with a power of 0.74 (x 2 /dof = 1.9) . 

5. CONCLUSIONS 

We see no evidence for the divergence predicted 
by the simple instanton model. Two possible con- 
clusions are reasonable. It is possible that the 
simulation is in some way insufficient to probe the 
physics from which this divergence derives it's ori- 
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2 /3 CJVt= 4 was taken to be 5.69 jij 



